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Introduction 


The program SCBUCKLE calculates the buckling loads and mode shapes of 
cylindrically curved, rectangular panels (see fig. 1). The panel is assumed to 
have no imperfections. SCBUCKLE is capable of analyzing specially 
orthotropic symmetric panels (i.e., A ]6 = A 26 = 0.0, Die - D 26 - 0.0 , Bjj = 0.0). 
The analysis includes first-order transverse shear theory, and is thus capable 
of modeling sandwich panels. The analysis supports two types of boundary 
conditions: either simply supported or clamped on all four edges. The panel 
can be subjected to linearly varying normal loads N x and N y in addition to a 
constant shear load N X y. The applied loads can be divided into two parts, a 
preload component, and a variable (eigenvalue-dependent) component. The 
analysis is based on the modified Donnell’s equations for shallow shells (see 
ref. 1). The governing equations are solved by Galerkin’s method. More 
details on the present analysis, including the limitations of shallow shell 
theory and important assumptions about the inplane boundary conditions, are 
provided in reference 2. 


Symbols 


a 

n 


An, A 12, A22, A 66 
A 16 , A26 
b 

Bij 

C x , Cy , C X y 


panel length (see fig. 1) 

amplitude coefficients of normal displacement series 
(see eqs. 8 and 9) 

orthotropic inplane stiffnesses of panel 
anisotropic inplane stiffnesses of panel 
panel width (see fig. 1) 

bending-stretching coupling matrix of panel ( i,j = 1, 2, 6) 
constant N x , N y , and N xy components, respectively 


/.0 ,»0 
Cjr, Cy f t X y 

Dq X > &Qy 

Du, D] 2 , D22 , D& 
Dw, D26 
@xz> @yz 



m, n 
rn 0 , n 0 

N x , Ny, N X y 

r 

tc 

tf 

x, y 
w 

X 


constant N x , N y , and N xy preload components, 
respectively 

transverse shear rigidities in the x- and y-directions, 
respectively (see eqs. 1-4) 

orthotropic bending stiffnesses of panel 

anisotropic bending stiffnesses of panel 

transverse shear stiffnesses in the x- and y-directions, 
respectively 

linear N x and N y components, respectively 
linear N x and N y preload components, respectively 

indices in assumed series solutions (see eqs. 8 and 9) 

number of terms in the x- and y-series, respectively 

longitudinal, transverse, and shearing stress 
resultants, respectively (see fig. 2) 

panel radius (see fig. 1) 

core thickness (see fig. 1) 

panel or facesheet thickness (see fig. 1) 

coordinate system (see fig. 1) 

buckling mode normal displacement (see eqs. 8 and 9) 
buckling eigenvalue (see eqs. 5, 6, and 7) 
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Note: Symbols in Courier font represent program variables. These variables 
are often a direct representation of physical or mathematical quantities (e.g., 
All represents An, Gxz represents Gxz, mO represents m 0 , etc.). 


Program Description 

The program is written in FORTRAN. The program consists of a main 
procedure and several subroutines, and was developed on a CONVEX C240 
computer under the UNIX operating system. The subroutine SYMGEP, 
which solves the symmetric eigenvalue problem, was obtained from the NASA 
Langley Research Center Computer Applications Branch. SYMGEP has its 
own documentation (see reference 3). Appendix A of this document includes 
the appropriate pages of reference 3. 

The program is capable of taking advantage of vectorizing compilers. By 
using the -02 option (i.e., local scalar optimization, global scalar 
optimization, and vectorization) in the CONVEX FORTRAN fc compiler, an 
almost tenfold reduction in CPU execution time was achieved. 


Input 

The program reads its input from the file inpt.dat. There are 31 lines in the 
file inpt.dat, with one entry per line. The program reads numerical data in 
free format. Dimensional variables may use any set of consistent units. 
However, all dimensional variables must use the same set of units. 

Appendix B contains an example set of input and output. 


A typical input file is: 


Test Case 42 
10.0 

! a 

- panel 

! title (32 characters maximum) 
length [Length] 

30.0 

! b 

- panel 

width [Length] 

60.0 

! r 

- panel 

radius [Length] 

0.32 

! tf 

- panel 

or facesheet thickness [Length] 

0.0 

! to 

- core i 

thickness [Length] 
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0 . 245888E+07 

! All 

- 

0 . 747557E+06 

! A12 

- 

0 . 245888E+07 

! A22 

- 

0.855662E+06 

! A66 

- 

0 . 209825E+05 

! Dll 

- 

0.637916E+04 

! D12 

- 

0.209825E+05 

! D22 

- 

0. 730165E+04 

! D66 

- 

0.5E+06 

! Gxz 

- 

0. 5E+06 

! Gyz 

- 

0.0 

! kxO 

- 

0.0 

! cxO 

- 

0.0 

! kyO 

- 

0,0 

! cyO 

- 

0.0 

! cxyO 

- 

6. 67E-02 

! kx 

- 

-1.0 

! cx 

- 

0.0 

! ky 

- 

0.0 

! cy 

— 

1.0 

! cxy 

- 

10 

! mO 

- 

10 

! nO 

- 

1 

! isc 

- 

l 

! iv 

- 


[Force/Length] 
[Force/Length] 
[Force/Length] 
[Force/Length] 
[Force -Length] 
[Force-Length] 
[Force-Length] 
[Force-Length] 


orthotropic inplane stiffness 
orthotropic inplane stiffness 
orthotropic inplane stiffness 
orthotropic inplane stiffness 
orthotropic bending stiffness 
orthotropic bending stiffness 
orthotropic bending stiffness 
orthotropic bending stiffness 
transverse shear stiffness [Force/Length A 2 ] 
transverse shear stiffness [Force/Length A 2] 
linear Nx preload component [Force/Length A 2] 
constant Nx preload component [Force/Length] 
linear Ny preload component [Force/Length A 2] 
constant Ny preload component [Force/Length] 
constant Nxy preload component [Force/Length] 
linear Nx component [Force/Length A 2 ] 
constant Nx component [Force/Length] 
linear Ny component [Force/Length A 2] 
constant Ny component [Force/Length] 
constant Nxy component [Force/Length] 
number of terms in x series (integer) 
number of terms in y series (integer) 
simply supported/clamped flag (integer) 
eigensolver flag (integer) 

- # of eigenvectors written to output file (integer) 


Comments to the right of the exclamation mark are not necessary for input; 
they are merely used for the purpose of documentation. The program will 
ignore the exclamation marks and comments when reading the input file. 

The quantities in square brackets indicate the dimensions of the variable. In 
the above example, forces are in units of pounds and lengths in units of inches. 
A brief explanation of the input variables follows: 


title - Title line for the input file. Any string of up to 32 characters in 
length is admissible. This input line must be present even if it is blank. 

a , b - Panel length and width, respectively. See figure 1. 

r - Panel radius. See figure 1. For a flat plate analysis r can be set to 
some very large number. 

t f - Panel or facesheet thickness. See figure 1. If the panel is of 

monolithic construction, t f is the total panel thickness. If the panel is 
of sandwich construction, t f is the thickness of a single facesheet; both 
facesheets are assumed to have the same thickness. Note that t f has 
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two interpretations. The variable tc (see next entry) determines the 
interpretation of tf. 

tc - Core thickness. See figure 1. If the panel is of monolithic 

construction, set tc = 0.0. Then tf is interpreted as the total panel 
thickness. If tc * 0.0, the panel is assumed to be of sandwich 
construction and t f is interpreted as the facesheet thickness (see 
previous entry). 

Note: The importance of t f and t c comes about in the calculation of 
the transverse shear rigidities Dq x and DQ y . If the panel is of 
monolithic construction, Dq x and Dq v are calculated by the formulas 

Dq x = | G xz tf (1) 

and 

DQy = ^ G yz tf (2) 

If the panel is of sandwich construction, Dq x and DQ y are calculated 
by the formulas 

D Qx = (3 ) 

and 

DQy = Gyz^-^t (4) 

J tc 

All, A12 , A22, A66 - Orthotropic inplane stiffnesses of panel. The 
panel is assumed to be specially orthotropic and thus, Aie = A 26 = 0.0. 

Dll, D12, D22, D66 - Orthotropic bending stiffnesses of panel. The 
panel is assumed to be specially orthotropic and thus, Djq = D 26 = 0.0. 
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Gxz , Gyz - Transverse shear stiffnesses in the x- and ^-directions, 
respectively. If the panel is of monolithic construction, these quantities 
are the transverse shear stiffnesses of the panel. If the panel is of 
sandwich construction, these quantities are the transverse shear 
stiffnesses of the core material. If no transverse shear deformation is 
desired in the analysis, Gxz and Gyz should be set equal to some large 
number. In this case setting Gxz and Gyz to be two or three orders of 
magnitude greater than their actual values is usually sufficient. 

kxO, cxO, kyO, cyO, cxyO, and kx, cx, ky, cy, cxy - Preload and 
load components. The stress resultants N x , Ny, and N X y define the 
loading on the panel. The longitudinal stress resultant N x is assumed 
to vary linearly with y and is given by: 

N x = k x y + c° x + X{k x y + c x ) (5) 

The transverse stress resultant N y is assumed to vary linearly with x 
and is given by: 


Ny = kyX + C*y + A [ky X + Cy) (6) 

The shear stress resultant N xy is assumed to be constant with respect 
to x and y and is given by: 

N X y — Cjcy + A C X y (7) 

fhe stress resultants N x and N y are positive in tension, while N xy is 
positive in its usual sense (see fig. 2). The program calculates the 
parameter A, which is the eigenvalue. Once A is known, the buckling 
stress resultants can be calculated from equations 5, 6, and 7. 

the variables kxO, cxO, kyO, cyO, cxyO, kx, cx, ky, cy, and cxy 

allow the user to define the loading on the panel. These variables are 
intended to be defined such that the loads on the panel have the desired 
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sign assuming a positive value of the eigenvalue, A. Variables with a 
zero (e.g., kxO, cxO, etc.) are the preload components. As can be seen 
from equations 5, 6, and 7, the loads defined by these preload 
components are not affected by the eigenvalue. The loads defined by the 
load components, kx, cx, ky, cy, and cxy, are affected by the 
eigenvalue. Although all of the preload components may be set equal to 
zero, at least one of the load components must be nonzero. 

The usage of these load components is illustrated by the following two 
examples. The first example is a uniform compressive N x load. The 
user sets cx - -1.0, and all other load and preload components equal 
to zero. The program will return a positive value of A as its lowest 
eigenvalue. The compressive buckling load, N x , can then be calculated 
from equation 5. This buckling N x will be negative. The second 
example is a N x load with a fixed gradient with respect to y but an 
unknown constant component. The user sets kxO equal to some 
nonzero value, cx — — 1 . 0, and all other load and preload components 
equal to zero. Again, once the program returns a value for A, the 
buckling load, N x , can be calculated from equation 5. 

The variables kx and ky allow the user to define linearly varying 
portions of the loads. However, the gradients of N x and N y with 
respect to y and x respectively may not be defined a priori if kx and/or 
ky are not equal to zero. These gradients also depend on A, which is a 
result of the program calculations, and thus not known in advance. If 
the gradients of N x and N y are to be specified a priori , it should be 

done through the variables kxO and kyO. 

Care should be exercised in choosing values for the preload components. 
If these components are sufficient by themselves to buckle the panel, the 
program will be unable to return a value for A. See the Output section 

for more information. 

m 0 , n0 - Number of terms in the assumed series in the x- and y- 
directions, respectively. Although mO and nO are two independent 
entries, the current implementation of the program requires them to 
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have the same value. Setting mO = nO = 10 will usually give 
acceptable results. The current implementation of the program 
requires mO and nO to have a value that is no greater than 20 due to the 
sizing of the arrays. These arrays can be resized as discussed in 
Appendix C. 

isc - Flag to specify boundary conditions. For simply supported boundary 
conditions set isc = X. For clamped boundary conditions set isc — 2. 

Note: Due to the limitations of shallow shell analysis, the present 
clamped boundary solution is accurate only for very small values of blr. 
As a result, this version of SCBUCKLE should not be used for curved 
panels with clamped boundary conditions. 

iv - Flag to specify whether or not to calculate the eigenvectors. If 

iv = 0, no eigenvectors are calculated. If iv = 1, all eigenvectors are 
calculated. 

nmodp • Number of eigenvectors to be written into the output file, starting 
with the eigenvector of the first eigenvalue. The program reads this 
line, but ignores it if iv = 0. Care should be used in selecting nmodp 
since large values result in voluminous output. The number of output 
lines generated by the program is greater than the product 
nmcdp*mQ*nO. 

Output 

The program writes its output to the file outp.dat. The program first echoes 
the input onto the output file to provide a mechanism for error checking. Next, 
all eigenvalues are written to the output file. There are mO * nO eigenvalues. 

All positive eigenvalues are printed in ascending order followed by the negative 
eigenvalues also in ascending order in their signed sense (i.e., not in their 
absolute value sense). If the loading is specified as described previously, the 
user will usually be interested in the first eigenvalue. If the program detects 
that any eigenvalue is less than the first one in the absolute value sense, a 
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warning message (including the mode and value of this lower eigenvalue) will 
appear before the eigenvalue listing. Finally, if iv = 1 and nmodp > 0, the 
eigenvectors of the first nmodp eigenvalues are written to the output file. 

These eigenvectors are the coefficients, a mn , of the normal displacement 
series. For simply supported boundary conditions this series is: 

mo no / \ 

w — X s (8) 


For clamped boundary conditions this series is: 


w 


V V [ |(wi • + 1) 7T X ) 

= L L a m n[cos [- ) - cos\ f 


m * 1 n = 1 


, ,(n - l)^y\ + l)*y 

• COS r — 1 - cos' 


(9) 


If the preload components are sufficient by themselves to buckle the panel, 
a warning message will be printed stating that one of the matrices is not 
positive definite. In this case no eigenvalues or eigenvectors are printed and 
execution is terminated. 


Program Execution 

The program is intended to be executed in the batch or background modes. 
Care should be taken in selecting mO and nO; making these quantities greater 
than 10 to 12 will greatly increase the execution time of the program. Because 
of the added complexity in the calculations, the solution for clamped boundary 
conditions will take from five to eight times longer than an equivalent solution 
for simply supported boundary conditions. 


9 


Appendix A: SYMGEP Documentation 


MATHEMATICAL AND STATISTICAL SOFTWARE AT LANGLEY— PART I 
HATE 03/01/87 REPLACES COPY OATEO 10/01/77 SECTION F2.7 

TITLE Subroutine SYMGEP PAGE 1 OF 7 


LANGUAGE: FORTRAN 

PURPOSE: To compute all the eigenvalues and eigenvectors of the 

symmetric generalized eigenvalue problem Ax * *8x where 
A is symmetric and B is symmetric, positive definite by 
the Cholesky-Wilkinson algorithm. 

USE: CALL SYMGEP(MAX,N,A,B,E,IV,WK,IERR) 

MAX An input integer specifying the first dimension 
of arrays A and B as stated in the dimension 
statement of the calling program. 

N An input integer specifying the order of A and B 

where 1 _<_ N <_ MAX. 

A An Input/output two-dimensional real array. 

Input A contains the symmetric input A matrix 
(only the full upper triangle of A need 
be supplied). 

Output A contains the eigenvectors normalized 
to unit length if IV - 1. If IV = 0, A 
is destroyed. The A array is 
dimensioned with variable dimensions in 
the subroutine. Therefore, A must be 
dimensioned in the calling program with 
first dimension MAX and second 
dimension at least N. 

B An input two-dimensional real array containing 

the symmetric, positive definite input B 
matrix. Only the full upper triangle need be 
supplied and only this part of the B array is 
not destroyed during the computations. The 
array B is dimensioned with variable dimensions 
in the subroutine. Therefore, B must be 
dimensioned in the calling, program with first 
dimension MAX and second dimension at least N. 
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MATHEMATICAL AMO STATISTICAL SOFTWARE AT LANGLEY— PART I 
OATE 03/01/07 REPLACES COPY DATED 10/01/77 SECTION F2.7 

TITLE Subroutine SYMGEP PAGE 2 OF 7 


OUTPUT 

INFORMATION: 


RESTRICTIONS: 

METHOD: 


E An output one-dimensional real array containing 

the. eigenvalues In ascending order. The array E 
Is dimensioned with variable dimensions In the 
subroutine. Therefore* E must be dimensioned in 
the calling program by at least N. 

IV An Input integer eigenvector option parameter. 

- 0 Compute all eigenvalues and no 
eigenvectors. 

= 1 Compute all eigenvalues and eigenvectors. 

WK An array for working storage. The array WK may 

be a mul tl -dimensional array of any type. The 
only requirement on WK is that it must occupy at 
least 2N locations. 

I ERR An output Integer error code 
= 0 Normal return. 

= j jtb eigenvalue has not been determined 

after 30 Iterations. 

= 7N+1 B is detected to be non-positive 
definite. 

Upon return, the calling program should test 
this parameter. 

The eigenvalues are stored in E by ascending order. 

The eigenvector associated with the i^ eigenvalue is 

found In the i ^ column of A. If the error return J is 
made, eigenvalues and eigenvectors are correct (but 
eigenvalues may be unordered) for indicts 1 ,2, . . . , IERR- 
1. If the error return 7M+1 is made, then no 
eigenvalues or eigenvectors are computed. 

None 

The Cholesky decomposition of B into LL T where L is 
lower triangular is first performed. The composition 

L"*AL"* is then performed, resulting in a symmetric 
matrix with the same eigenvalues as the original matrix 
system. The eigenvalues and eigenvectors of this 
symmetric matrix are found by the explicit QL 
algorithm. These eigenvectors are then transformed into 
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MATHEMATICAL AND STATISTICAL SOFTWARE AT LANGLEY— PART I 


DATE 03/01/87 REPLACES COPY DATED 10/01/77 SECTION F2.7 

TITLE Subroutine SYMGEP PAGE 3 OF 7 


ACCURACY: 


REFERENCES: 


the eigenvectors of the original matrix system. The 
eigenvectors are normalized such that x T Bx = 1. 

This subroutine Is numerically stable; 1.e. t each 
computed eigenvalue and Its corresponding eigenvector 
are exact for a matrix problem very close to the 
original problem. The accuracy of the eigenvalues Is 
dependent upon the condition of 8 with respect to 
Inversion. If B Is well -conditioned, then the largest 
eigenvalue Is computed to approximately 12 significant 
figures and the smaller eigenvalues will suffer from 
absolute errors which are no larger. The accuracy of 
the eigenvectors Is dependent on their Inherent 
sensltlvl ty. 

1. Martin, R. $.; and Wilkinson, J. H.: Reduction of 

the Symmetric Elgenproblem, Ax = XBx and Related 
Problems to Standard Form. Mumer. Math., Bd. 11, 
I960, pp. 99-110. 

2. Martin, R. S.; Relnsch, C.; and Wilkinson, J. H.: 
Householder's Trldlagonallzatlon of a Symmetric 
Matrix. Mumer. Math. Bd. 11, 1968, pp. 181-195. 

3. Bowdler, M.; Martin, R.S.; Relnsch, C.; and 

Wilkinson, «J. II.: The QR and Ql Algorithms for 

Symmetric Matrices. Mumer. Math., Bd. 11, 1968, 
pp. 293-306. 


All the references are reprinted In Handbook for 
Automatic Computations, Volume II, Linear Algebra by J. 
II. Wilkinson and C. Relnsch, Sprlnger-Verlag, 1971. 


STORAGE : 

2623 0 


SUBPROGRAMS 

USED: 

SYM()I. (r?.l) 

73, 


QX2138 

426, 


QXZ139 

351, 


QXZ140 

304, 


QXZ155 

277, 


QXZ156 

110, 
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HATMEMATICAL WIO STATISTICAL SOFTWARE AT LANGLEY-PART I 


DATE 03/01/07 REPLACES COPY OATED 04/30/80 
TIILE Subroutine SYMGF.P 

QXZ168 264 0 
QXZ262 17 8 
QXZ263 31g 


SECTION E2.7 
PAGE 4 OF 7 


FORTRAN ABS. AMAX1, AMIN1 , IABS, SIGN, SQRT 

FUNCTIONS: 


OTHER CODING 
INFORMATION: 


If B Is Ill-conditioned with respect to Inversion 
subroutine CSQZ (F2.8) should be used for more accurate 
results. 


EXAMPLE 


If a user wanted to compute all the eigenvalues and 
eigenvectors of a 5 by 5 symmetric generalized 
elgenproblem, then the user must store at least the full 
upper trlahgle of the A and B matrices In two arrays, 
let A be stored In an array denoted A and dimensioned 5 
by 5 and let B be stored In an array denoted B and 
dimensioned 5 by 5. (Note that A and B must have the 
same first dimension). Then the elements of the A 
matrix, denoted a^j, and of the B matrix, denoted b f j, 

are stored as follows: 


A 


1 B 

i 

- 

a ll a 12 a l3 a 14 a i5 


b ll b 12 b 13 b 14 b 15 

x a 22 a 23 a 24 a 25 


x b 22 b 23 b 24 b 25 

x x a 33 a 34 a 35 


x x b 33 b 34 b 35 

X x x a 44 a 45 


x x x b 44 b 45 

X x x * a 55 


_x x x x b 55 


where x denotes elements which are Ignored by the 
subroutine. Thus, the full A and B matrices may be 
stored In A and B respectively, but only the full upper 
triangle will be used as Input. Let the array In which 
the computed eigenvalues are to be stored be denoted by 
E and dimensioned by 5. (E must be dimensioned by at 
least 5). Also, suppose that there exists an array 
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NAlflEJVU I CAL AND STATISTICAL SOFTWARE AT LANGLEY— PART I 


DATE 03/01/87 REPLACES COPY DATEO 01/01/85 SECTION F2.7 

TITLE Subroutine SYMGEP PAGE 5 OF 7 


denoted W which Is dimensioned 3 by 2 by 2 and whose 
contents may be destroyed. Since W occupies 12 
locations and 10 locations are needed for a work array, 
W may be used for this array. Then, a FORTRAN program 
performing this task could be as follows: 

PROGRAM TTTSYMG( OUTPUT , TAPE6=0UTPUT ) 

THIS ROUTINE COMPUTES THE COMPLETE EIGENVECTOR 
SYSTEM OF THE SYMMETRIC GENERALIZED EIGENVALUE 
PROBLEM BY THE CHOLESKY-WILKINSON ALGORITHM. 

DIMENSION A(5,5),E(5),B(5,5),W(3,2,2) 

DATA ( (A( I ,J) , 1*1 ,5) ,J*1 ,5)/lO. ,2. ,3. , 1 . , 1 . ,2. , 

■r I2.,l.,2.,l.,3.,l.,ll.,l.,-l.,l.,2.,l.,9.,l., 

+ t.,1.,-1.,1.,15./ 

DATA ( (B( I ,J) , 1*1 ,5) ,.J*1 ,5)/12. ,1 . ,-l . ,2 . ,1. , 1 . , 

♦ l 4 .,l. l -l.,l.,-l.,t., 16 .,-l.,l., 2 .,-l.,-l., 12 .,-l., 

♦ l.,l.,l.,-l.,U./ 

MAX * 5 

N = 5 
IV * 1 

WRITE (6,896) 

WRITE (6,899) ( (A( I ,J) ,J*1,5) , 1 = 1,5) 

WRITE (6,898) 

WRITE (6,899) ( (B( I ,J) ,J*1 ,5) , 1*1 ,5) 

CALL SYMGEP(MAX,N,A,B,E,IV,W,IERR) 

IF (IERR.EQ.O)GO TO 100 
WRITE (6,900) IERR 

896 FORMAT!/, 1X.37HTHE ORIGINAL SYMMETRIC REAL MATRIX A:) 

898 FORMAT!/, IX, 37HTIIE ORIGINAL SYMMETRIC REAL MATRIX B : ) 

899 FORMAT!/, 5(5X,F5.1)) 

STOP 

100 CONTINUE 

WRITE (6,903) 

WRITE (6,904) (E(K),K*1,5) 

WRITE (6,905) 

WRITE (6,906) ( (A( I ,J) , 1 = 1 ,5) ,J = 1 ,5) 

WRITE (6,900) IERR 

900 FORMAT (/, IX, 1 7HTIIE ERROR CODE IS, IX, 12) 

903 FORMAT ( / , l X , 20HTIIE EIGENVALUES ARE:) 

904 FORMAT!/, 1X.E20. 13) 

905 FORMAT!/, IX, 20HTIIE NORMALIZED EIGENVECTORS:) 

906 FORMAT (/,1X,5(1X,E20.13,/)) 

STOP 

ENO 
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MATHEMATICAL AMO STATISTICAL SOFTWARE AT LANGLEY— PART I 
DATE 03/01/87 REPLACES COPY DATED 10/01/77 SECTION E2.7 

TITLE Subroutine SYMGEP PAGE 6 OF 7 


If the j ^ eigenvalue Is denoted by ej and vjj denotes 

the 1^^ component of the eigenvector associated with the 

j th eigenvalue, then the E and A arrays contain the 
following elements upon a normal return: 


e l 


v ll 

v 12 

v 13 

v 14 

v 15 

p ? 


v 2 1 

v 22 

v 23 

v 24 

v 25 

*3 

A - 

v 3 l 

v 32 

v 33 

v 34 

v 35 

p 4 


V 41 

v 42 

v 43 

v 44 

v 45 

p 5 


V 51 

v 52 

v 53 

v 54 

< 

cn 

cn 

t - 


This Is shown by the following output from program 
TTTSYMG 

THE ORIGINAL SYMMETRIC REAL MATRIX A: 


10.0 

2.0 

3.0 

1.0 

1.0 

2.0 

12.0 

1.0 

2.0 

1.0 

3.0 

1.0 

11.0 

1.0 

- 1.0 

1.0 

2.0 

1 .0 

9.0 

1.0 

1.0 

1.0 

- 1.0 

t.O 

15.0 

ORIGINAL 

SYMMETRIC 

REAL MATRIX 

B : 


12.0 

1.0 

- 1 .0 

2.0 

1.0 

1.0 

14.0 

1.0 

- 1 .0 

1.0 

- 1.0 

1.0 

16.0 

- 1.0 

1.0 

2.0 

- 1.0 

- 1.0 

12.0 

- 1.0 

1.0 

1.0 

1.0 

- 1.0 

11.0 


THE EIGENVALUES ARE: 
.4327872110170E400 
.6636627483923E400 
. 9438590046684E ♦ 00 
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MATHEMATICAL AMO STATISTICAL SOFTWARE AT LANGLEY— PART I 


DATE 03/01/87 REPLACES COPY DATED 10/01/77 SECTIOM F2.7 

TITLE Subroutine SYMGEP PAGE 7 OF 7 


. 110928454001 7E+01 

. 1192353232513E »01 

THE MORHAI.IZEI) EIGENVECTORS: 

- 54 7 1 53885036 IE *^00 
219182 72 16828E +00 
-.61 19216001 273E+00 
, 11501 35612106E +00 
168601 1052691E*00 

.2972211326819E *-00 
.5189568898292E«00 
-.1251321 101 336E *00 
-.6552891292303E+00 
. 1276690116057E-01 

-.6238090551 130E^0 
.51731161081 72E *00 
-.2135203290675E+00 
.11 73 1286252 5 IE *00 
28952691 98812E +00 

.5333052571006E+00 
. 5318371066823E >00 
. 1 5138897 18681E * 00 
.17111 31667187E+00 
. 288069651 7888E - 0 1 

-. 2391 72806351 3E* 00 
. 5360233565381E -01 
- . 20R9937U21252E +00 
.26976101 13650E *00 
. 90 7 375112608 IF. +00 

THE ERROR CODE IS 0 


SOURCE: NASA, LaRC, Computer Applications Branch. 

QUESTIONS ON THE USE OF TIMS PROGRAM SHOULD BE DIRECTED TO THE ACD USER 
CONSULTATION OFFICE, EXT. 3518. 
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Appendix B: Example Input and. Output 


inptdat 


Test Case 42 
10.0 

! a 

- 

30.0 

! b 

- 

60.0 

! r 

- 

0.32 

! tf 

- 

o 

o 

! tc 

- 

0 . 245888E+07 

! All 

- 

0.747557E+06 

! A12 

- 

0 . 245888E+07 

! A22 

- 

0 . 855662E+06 

! A66 

- 

0 . 209825E+05 

! Dll 

- 

0 . 637916E+04 

! D12 

- 

0 . 209825E+05 

! D22 

- 

0 . 730165E+04 

! D66 

- 

0.5E+06 

! Gxz 

- 

0.5E+06 

! Gyz 

- 

0.0 

! kxO 

- 

0.0 

! cxO 

- 

0.0 

! kyO 

- 

0,0 

! cyO 

- 

0.0 

! cxyO 

- 

6. 67E-02 

! kx 

- 

-1.0 

! cx 

- 

0.0 

! ky 

- 

0.0 

! cy 

- 

1.0 

! cxy 

- 

10 

! mO 

- 

10 

! n0 

- 

1 

! isc 

- 

1 

! iv 

- 

1 

! nmodp 


! title (32 characters maximum) 
panel length [Length] 
panel width [Length] 
panel radius [Length] 
panel or facesheet thickness [Length] 
core thickness [Length] 

orthotropic inplane stiffness [Force/Length] 
orthotropic inplane stiffness [Force/Length] 
orthotropic inplane stiffness [Force/Length] 
orthotropic inplane stiffness [Force/Length] 
orthotropic bending stiffness [Force-Length] 
orthotropic bending stiffness [Force-Length] 
orthotropic bending stiffness [Force-Length] 
orthotropic bending stiffness [Force-Length] 
transverse shear stiffness [Force/Length A 2 ] 
transverse shear stiffness ( Force/Length A 2 ] 
linear Nx preload component [Force/Length A 2 ) 
constant Nx preload component [Force/Length] 
linear Ny preload component [Force/Length A 2 ] 
constant Ny preload component [Force/Length) 
constant Nxy preload component (Force/Length] 
linear Nx component [Force/Length A 2] 
constant Nx component [Force/Length] 
linear Ny component [Force/Length A 2] 
constant Ny component [Force/Length] 
constant Nxy component [Force/Length] 
number of terms in x series (integer) 
number of terms in y series (integer) 
simply supported/clamped flag (integer) 
eigensolver flag (integer) 

- # of eigenvectors written to output file (integer) 


outadat 


+f +f ++++++ 

+SCBUCKLE+ 

+ + + + + 

Case: Test Case 42 

*** Panel Definition *** 

a - 0 . 100000E+02 [Length] 

b - 0 . 300000E+02 [Length] 

r - 0 . 600000E+02 (Length) 

tf - 0 . 320000E+00 [Length] 

tc - 0 . OOOOOOE+OO [Length] 

*** Load Definition *** 

kxO - 0 . 000000E+00 [Force/Length A 2] 

cxO - 0 . OOOOOOE+OO (Force/Length] 

kyO - 0 . 000QQ0E+00 [ Force/Length A 2 ] 

cyO - 0 ♦ OOOOOOE+OO [Force/Length] 

cxyO - 0. 000000E+00 [Force/Length] 

kx - 0 . 667000E-01 [ Force/Length A 2 ] 
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cx - -0 . 100000E+01 [Force/Length] 
ky * 0 . 000000E+00 [ Force/Length A 2 ] 

cy - 0 , 000000E+00 {Force/Length] 

cxy * 0 ♦ 100000E+01 [Force/Length] 

»** Inplane Constants *** 

All - 0 . 2 45888E+07 [Force/Length] 

A12 - 0 . 747557E+06 ( Force/Length J 

A 22 « 0 . 2 45888E+07 {Force/Length] 

A66 - 0 . 855662E+06 {Force/Length] 

*** Bending Constants *** 

DU - 0.209825E+05 (Force-Length 1 

D12 • 0 . 63791 6E+04 [Force-Length J 

D22 - 0 .2098252+05 (Force-Length] 

D66 - 0 . 730165E+04 [Force-Length] 

*** transverse Shear Constants *** 

Gxz - 0.500000E+06 [Force/Length A 2 ] 

Gyz - 0 . 500000E+06 (Force/Length A 2 ] 

*** Program Execution and Output *** 
mO - 10 

nO - 10 

isc - 1 

iv - 1 

nmodp - 1 


+ + M- + + +*+++++*+ 

+ + + Results +++ 

+ *■ + + + + + + + + + + f + + 

Warning 

The first eigenvalue does not have the lowest absolute value 

Mode 100, lambda * -0 . 907929E+04 


★** Eigenvalues *** 


Mode lambda 

1 0 . 908700E+04 

2 0 . 103833E+05 

3 0 . 1542 61E+05 

4 0 . 159696E+05 

5 0.207969E+05 

6 0 . 224853E+05 

7 (5 . 23382BE+05 

8 0.26226 2E+Q5 

9 0 . 2 98956E+05 

10 0 . 310562E+05 

11 0 . 3281 90E+05 

12 0 . 356992E+05 

13 0 . 402 612E+05 

14 0 . 4 4172 4 E + 05 

15 0 . 4 602 30E+G5 

16 0 . 483892E+05 

17 Q. 541446E*-05 

18 0 . 589957E+05 

19 0 . 623080E+05 

20 0 . 6320 10E+05 

21 0 . 674 64 8E+05 
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22 0 . 690096E+05 

23 0 . 733284E+05 

24 0 . 757669E+05 

25 0 . 789729E+05 

26 0 . 825614E+05 

27 0 . 8617 92E+05 

28 0 . 883095E+05 

29 0 . 91301 1E+05 

30 0 . 934 032E + 05 

31 0 . 1G757 9E + 06 

32 0 . 1 15188E+06 

33 0 . 126603E+06 

34 0 . 1 3088 1 E+06 

35 0 . 14 4 7 69E+06 

36 0 . 151 668E+06 

37 0.16144 6E+06 

38 0 . 168785E+06 

39 0 . 200883E+06 

40 0 . 215595E+06 

41 0 . 226502E+06 

42 0 .2282 92E+06 

43 0 . 2 75 350E+0 6 

44 0 » 284 34 1E+06 

45 0 . 4 48697E + 06 

46 0 . 4 67881E+06 

47 0 , 65972 5E +06 

48 0 . 698243E + 06 

49 0 . 817984E +06 

50 0 . 86128 1E+06 

51 -0.852285E+06 

52 -0 . 809542E+06 

53 -0 . 691 4 39E+06 

54 -0 . 653236E+06 

55 -0.463890E+06 

56 -0.444889E+06 

57 -0 . 283333E+06 

58 -0 . 274285E+06 

59 -0.227505E+06 

60 -0 . 2255 92 E+06 

61 ~0 . 2 1427 3E +06 

62 -0 . 199933E+06 

63 -0 . 1684 27E + 06 

64 -0 . 1 61072E+0 6 

65 -0 . 151252E *-06 

66 -0 . 144397E+06 

67 -0 . 13063 1E+06 

68 -0 . 126358E+06 

69 -0.115011E+06 

70 -0 . 1074 18E+06 

71 -0 . 932592E+05 

72 -0 . 911722E+05 

73 -0 . 881 7 34 E +05 

74 -0.860516E+05 

75 -0 . 824 105E-+05 

76 -0 . 788320E+05 

77 -0 . 756697E+05 

78 -0 . 7313 93E+05 

79 -0 . 68874 6E+05 

80 -0.673460E+05 

81 -0 . 631025E+05 

82 -0 . 622 398E+05 

83 -0 . 589280E+05 

84 -0.540839E+05 

85 -0 . 4 83303E+05 
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86 

-0 . 4 59742E+05 

87 

-0 . 4 4 1282E+05 

88 

-0 , 4 02225E + 05 

89 

-0 . 356657E+05 

90 

-0 * 327877E+05 

91 

-0 . 310288E+05 

92 

-0 . 298654E+05 

93 

-0 . 262048E+05 

94 

-0 . 233619E+05 

95 

-0.224 666E+05 

96 

-0 . 207773E*05 

97 

-0 . 159561E+05 

98 

-0454127E+05 

99 

-0403753E+05 

100 

-0 . 907 928E+04 


*** Eigenvectors *** 


Mode 1 


m n Amn [Length] 


1 

1 

1 

1 

l 

1 

1 

l 

1 

1 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 
A 
A 
A 
A 
A 
A 
A 
A 
A 

5 


X -0 . 28660 1E+00 

2 -0 . 6522 03E + 00 

3 -0.6H955E+00 

4 -0. 187865E+00 

5 0.24 98 68E-0 1 

6 0.36X405E-0X 

7 0 . 201525E-01 

8 0404224E-01 

9 0 . 591865E-02 

10 0.363X56E-02 

X 04Q8739E+00 

2 0 ♦ 057 X X0E-0X 

3 -Q , 94 6775E-01 

4 -0478772E + 00 

5 -0 . 108622E+00 

6 -0.446994E-01 

7 -Q.179662E-0X 

8 -0.84 7 606E-02 

9 -0.447371E-02 

10 -0.224728E-02 

X 0.441720E-02 

2 0 , 233624E-01 

3 0.327385E-01 

4 0 . 993520E-02 

5 -G.1QQ144E-01 

§ - 04299385*01 

7 -0 . 926790E-02 

8 -Q.5939X0E-02 

9 -0 . 394 7 69E-02 

10 -0 . 255752E-02 

1 0 . 254 005E-02 

2 0 , 21 1561E-02 

3 -0 , 2226B4E-02 

4 -0 . 401 1 06E-G2 

3 -0 . 247 707E-02 

6 -0466062E-02 

7 -0 . 135056E-02 

8 -0 . 117728E-02 

9 -0 . 935713E-03 

10 -0 . 789628E-03 

1 0432101E-03 



5 

2 

0 . 15082 7E-02 

5 

3 

0 . 24 322 9E-02 

5 

4 

0 . 728031E-03 

5 

5 

-0 . 10005 7E-02 

5 

6 

-0 ♦ 124696E-02 

5 

7 

-0 . 101230E-02 

5 

8 

-0 ♦ 768483E-03 

5 

9 

-0 , 602899E-03 

5 

10 

-0.485602E-03 

6 

1 

0.417151E-03 

6 

2 

0.377364E-03 

6 

3 

-0 . 384930E-03 

6 

4 

-0.751984E-03 

6 

5 

-0.477391E-03 

6 

6 

-0 . 317430E-03 

6 

7 

-0.261496E-03 

6 

8 

-0.247351E-03 

6 

9 

-0 , 21 3723E-03 

6 

10 

-0 . 207705E-03 

7 

1 

0.184008E-04 

7 

2 

0 . 316782E-03 

7 

3 

0 . 535752E-0 3 

7 

4 

0 . 160534E-03 

7 

5 

-0 . 236204E-03 

7 

6 

-0 . 298661E-03 

7 

7 

-0 . 250287E-03 

7 

8 

-0.199534E-03 

7 

9 

-0 . 1 64 887E-03 

7 

10 

-0. 144380E-03 

8 

1 

0. 126053E-03 

8 

2 

0. 117965E-03 

8 

3 

-0 . 118395E-03 

8 

4 

-0.23 92 60E-03 

8 

5 

-0 . 1534 B7E-03 

8 

6 

-0 . 101536E-03 

8 

7 

-0 . 8351 94E-04 

8 

8 

-0 . 81 174 6E-04 

8 

9 

-0 . 720956E-04 

8 

10 

-0 . 742 67 4E-04 

9 

1 

0 . 503171 E— 0 5 

9 

2 

0 . 109984E-03 

9 

3 

0.189801E-03 

9 

4 

0 . 569054E-04 

9 

5 

-0 . 86037 5E-04 

9 

6 

-0 . 10937 3E-03 

9 

7 

-0.92 7 95 IE-04 

9 

8 

-0 . 754565E-04 

9 

9 

-0.637054E-04 

9 

10 

-0 . 57 94 67E-04 

10 

1 

0 . 525962E-04 

10 

2 

0 . 4984 95E-04 

10 

3 

-0.4947 67E-04 

10 

4 

-0 . 100961E-03 

10 

5 

-0 . 64 6068E-04 

10 

6 

-0 . 428839E-04 

10 

7 

-0 . 355222E-04 

10 

8 

-0 . 35317 9E-04 

10 

9 

-0 . 318578E-04 

10 

10 

-0 . 34 3392E-04 
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Appendix C: Array Resizing 


The user may wish to resize the arrays in the program for use on smaller 
or larger computers. The arrays in the main program can be resized by 
changing the values of mmax, nmax, and mnmax. Although mnmax and nmax 
are independent quantities, they must have the same numerical value. The 
quantity mnmax must be set equal to the product mmax* nmax. If these 
quantities are changed, the arrays in the subroutines SSBUCK and CCBUCK 
must also be resized to agree with the new values of mmax, nmax, and mnmax. 
The arrays in these subroutines must be resized as follows: 


Subroutine SSBUCK array re sizing 


Lmn 

PI through P13 
Ql, Q2, and Q4 
dl and gl 
d2 and g2 
bb 


-> 

-> 

-> 

-> 

-> 

-> 


(nmax, nmax) 
(mmax, nmax) 
(mmax, nmax) 
(mmax, nmax) 
(nmax, nmax) 
(mnmax, mnmax) 


wk -> (2, mnmax) 

hroutine CCBUCK array 

resizing 

R1 through R4 

-> 

(mmax, nmax) 

SI through $4 

-> 

(mmax, nmax) 

P 

-> 

(48, nmax, nmax) 

Ql through Q4 

-> 

(nmax, nmax) 

Q9 through Q12 

-> 

(nmax, nmax) 

dl through d4 

-> 

(rrmax+2 , nmax+2 ) 

VI through VI 0 

-> 

(mnmax, mnmax) 

tempi through tempi 9 

-> 

(mnmax, mnmax) 

bb 

-> 

(mnmax, mnmax) 

wk 

-> 

(2, mnmax) 


where the numeric values of mmax, nmax, and mnmax are used. 
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Figure 1. Panel Geometry. 
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Figure 2. Stress Resultant Sign Convention, 
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